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We present a simple and direct proof that the exchange-correlation hole, and therefore the 
exchange-correlation energy, in a polarized insulator is not determined by the bulk density alone. It 
is uniquely characterized by the density and the macroscopic electric polarization of the dielectric 
medium. 
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In a famous paper in 1964 Hohenberg and Kohn (!]] 
established the basic principle of density functional the- 
ory (DFT): any property of an interacting time-reversal 
invariant many-particle system is uniquely determined 
as a functional of its ground state density no(r). The 
proof rests upon the demonstration that a given no(r) is 
consistent with only one external scalar potential v ex t(r) 
(except for an arbitrary constant), which determines the 
entire Hamiltonian since kinetic and particle-interaction 
{V e - e ) terms are known from the fundamental constants 
h, m e , and e. Together with the Kohn-Sham (KS) || 
ansatz that the interacting ground state density can be 
represented by non-interacting fermions in an effective lo- 
cal potential, this has provided the foundation for much 
of the current theoretical calculations in atomic, nuclear 
and condensed matter physics || . 

In the past DFT has often been applied uncritically to 
condensed matter, assuming implicitly that the proper- 
ties of bulk matter in the thermodynamic limit are func- 
tional only of the density in the bulk. However, the 
Hohenberg-Kohn (HK) proofs are only for the entire den- 
sity (including the surface) and they provide no guaran- 
tee that the functionals of the density are well-defined in 
the thermodynamic limit. In particular, the use of Born- 
von Karman (BvK) boundary conditions may introduce 
additional considerations not included in their original 
analysis. Although it is well-known that the long-range 
Coulomb terms must be treated carefully and are not 
absolutely convergent in the thermodynamic limit, it has 
been implicitly assumed that the key quantity describing 
many-body effects in the KS approach, the exchange- 
correlation (x-c) energy E xc , is an intrinsic functional of 
the bulk density. Recently, however, Gonze, Ghosez, and 
Godby (GGG) g) have argued that in an insulator E xc 
must also be a function of the macroscopic (average) elec- 
tric polarization in the limit of an infinite system. This 
has led to much controversy and proposals concerning 
the nature of exchange and correlation and the KS 
functional [plftl in insulators. 



is modified in a manner not described solely by the den- 
sity in the bulk region of the crystal, even though the 
hole around an electron in the bulk is localized to the 
bulk region. We show that the hole is determined only 
if the bulk macroscopic electric polarization is specified 
in addition to the bulk density. This can be interpreted 
as the "polarization" of the x-c hole HQ, providing the 
physical mechanism leading to the generalized form of 
the KS equations required to have an in-principle-exact 
description of polarized extended matter |^|,(?J . 

In a dielectric with a macroscopic electric field E mac , 
the energy can be expressed in terms of the polarization 
|iof . In the case of a macroscopic region of a periodic 
crystal of volume £7, external effects can be incorporated 
by terms of the form (fi = m e = — e = 1 in a.u.) 



E, 



dr v ext (r) n(r) 



(1) 



where the field couples to the macroscopic polarization 
Pmac i and the second term is explicitly restricted to only 
the periodic part of the external potential v ex t(r) cou- 
pled to the electronic (number) density n(r), which is 
periodic. The key point for our purposes is that P mac is 
not directly determined by the density in Q. Although 
variations in the polarization density field are related to 
n(r) by V • P(r) = +n(r), its average value P mac can 
be found only as a change from a reference polarization 
P^, ac , in terms of an integrated polarization current that 
flows through the volume f2 fTi"|| . This result for P mac is a 
consequence of the BvK boundary conditions used to de- 
scribe the extended system Jl2| ; it is not needed in finite 
systems and was not considered in the original analysis of 
HK who, implicitly, assumed open (vanishing) boundary 
conditions. Thus, the bulk density n(r) and changes in 
polarization <5P mac — P mac — P° lac are independent ex- 
perimentally measurable quantities. Furthermore, each is 
uniquely determined by the bulk ground state wave func- 
tion, with <5P mac given by its Berry's geometric phase 



In this Letter we provide a simple direct demonstra- 
tion that the x-c hole in a polarized insulating crystal 



The case of zero macroscopic field, E mac = 0, is of spe- 
cial interest because it can be considered the standard 
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reference state of the insulator: then its total energy is 
not explicitly dependent upon P mac 0, and the original 
analysis of HK is sufficient to conclude that all bulk prop- 
erties must be Junctionals only of the bulk density n(r); 
in particular, this includes the ground state spontaneous 
polarization P^, ac . On the other hand, if E mac is not 
zero, the same reasoning requires that all physical quan- 
tities must be functions of SP nmc as well as functionals 
of the bulk density. Therefore, the intrinsic description 
of a dielectric material in the presence of an electric field 
requires generalization of the HK theorems [Q . 

In the KS approach, E xc is defined to be the difference 
between the internal energy of the interacting many-body 
system and that of a non-interacting system of electrons 
of the same density n(r). The generalization for a po- 
larized system requires that the non-interacting system 
also has the same zero-field polarization P^, ac and change 
<5Pmac in the presence of a field. As shown in Ref. 0, 
E xc is related to the interacting N-particle wave function 
^ by a generalized constrained Levy formulation: 

1 N 

E xc [n; P mac ] = min (*|--VvJ + V e - e \9) 

f^n r ,P mac 2 f— ' 

i—l 

-T ind [n;P nmc ] - E H [n] , (2) 

where Ti ne i is the kinetic energy of independent fermions 
and Ejj is the Coulomb Hartree energy JtJ. This may 
also be expressed in terms of the Coulomb interaction 
between the electron density at point r and the density 
at r + uof the x-c hole n xc (r, r + u) defined by a coupling 
constant integration [|| over all values of the interaction 
A from A = to its full strength A = 1, 

E xc = \ [ dr n(r) / du -^-n xc (r, r + u) 
z Jn J o l u l 

N f 1 

= 2 y a du O Ke( U )) > ( 3 ) 

where (n xc (u)) is the density- weighted average of the 
x-c hole. Equivalently, E xc can be written in terms of 
the coupling constant averaged (symmetric) pair corre- 
lation function g x through the relation n xc (r, r + u) = 
Jo dX [g x {r, r + u) - l] n(r + u). 

The energy E XC /~N can be interpreted as the electro- 
static interaction between an electron of charge —1 and 
the charge distribution — (n xc (u)), the total charge of 
which is +1, as required by a sum rule. Thus E xc is de- 
termined by the shape of (n xc (u)). The shape is also re- 
stricted by other considerations, e.g., it is straightforward 
to show that the density- weighted average of n xc (r, r + u) 
over the entire volume of any finite system f2 must be 
symmetric in u, i.e., (n xc (u)) = {n xc (— u)}. This fol- 
lows from the fact that correlations among electrons are 
symmetric in particle exchanges. In addition, it can be 
shown that averages over r in a cell of a periodic region 
of a material also lead to a symmetric (n xc (u)) as long as 



the hole is localized in the relative coordinate u, so that 
r + u is also limited to the periodic region. This applies 
to a crystal in an electric field, where the (metastable) 
ground state wave function is invariant to lattice trans- 
lations even though the Hamiltonian is not. In this case, 
however, we note that other averages of the x-c hole, such 
as the simple average over position, J Q dr n xc (r,r + u), 
are not symmetric. 

Now let us consider the case of an insulator in an elec- 
tric field (or, as we will see, the case of a spontaneously 
polarized insulator) . The goal is to establish that the x-c 
hole is both a functional of the bulk density and a func- 
tion of the polarization; for this it is sufficient to consider 
only the non-interacting case where we can solve exactly 
for the exchange (x) hole. It then follows from the fact 
that E xc can be expressed as the coupling constant av- 
eraged form, for which the non-interacting x hole is one 
limit, that the same conclusions must apply to the exact 
x-c hole. We consider a one-dimensional (Id) insulat- 
ing crystal with potential energy V(x) = Vq(x) + Vl(x), 
where Vq(x) is periodic in the crystal lattice constant 
a, and Vl(x) is an additional long wavelength poten- 
tial with periodicity L = Ma, M > 1. In order to 
mimic a constant electric field near x — we have cho- 
sen V L (x) = V b [|§ sin(i) - ^ sin(3x) + gfg sin(5i)] , 
with x = 2ttx/L, which is a symmetric form that is 
oc x + 0(x 7 ). Figure [l] shows the potential energy V(x) 
and the change in density, An(x) — n(x) — n(x)v b =o, in 
the region around x — 0. Although the potential clearly 
has a linear component, the density is very nearly peri- 
odic. In the middle panel is also shown the change in the 
periodic density generated by a potential V^ lt (x) which 
has periodicity a chosen to reproduce the density of the 
central part of the supercell. Since the curves are essen- 
tially indistinguishable, this illustrates the point made by 
GGG that the same periodic density can be generated by 
two potentials differing by more than a constant. 

Consider the density-weighted average x hole around 
an electron in the periodic region around x — of the 
supercell. We define this hole to be (n^^u)) and the 
change from the zero applied field hole to be A(ni 1 ' ) (?i)), 
which is readily expressed in terms of the eigenstates of 
the supercell. We have verified that A{n x (u)) does not 
depend upon the unit cell chosen for the average in a 
region around the origin where the electric field is effec- 
tively constant. Also, the lowest order change in the hole 
must be quadratic in the field strength since the zero- 
field system has a center of inversion symmetry, and we 
have plotted in Fig. [l] the change per unit field squared. 
The result is compared with the corresponding change for 
the system with the same density and potential V' lt (x), 
which is defined to be A{n x (u)). It is apparent from the 
figure that A(ni 1 ' ) (u)) is more extended than A(n x 2 \u)) . 
Since the two holes are different, we have established that 
the density is not sufficient to determine the x hole, and 
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using the reasoning given above we conclude that E xc is 
not determined solely by the bulk periodic density. 

It is straightforward to prove that the hole is deter- 
mined by the bulk density and the macroscopic polariza- 
tion. This follows from the same reasoning as the original 
HK arguments except that now the external terms in the 
energy (Eq. d3)) involve both the bulk n(r) and P mac - In 
our present example, we have calculated the polarization 
in the two cases to show explicitly that they are differ- 
ent. The average polarization P^lc = 3.44 x 10 -3 in 
cells near the origin of the supercell can be found from 
the density in the supercell using dP(x)/dx = +n(x), 
whereas the polarization Pm ac = 2.11 x 1CP 3 of the sys- 
tem with the fitted periodic potential can be found from 
the Berry's phase expressions |l3|,|l4|]. To our knowledge, 
this is the first explicit demonstration of the polarization 
dependence of exchange, and therefore of the x-c energy. 

It is instructive to consider the Clausius-Mossotti limit 
of localized, non-overlapping wave functions in each cell. 
Then the two different ways of calculating P mac lead to 
the same answer: the dipole moment per unit length of 
each isolated unit. The x hole is also localized in one 
cell, so that all properties are determined by the density 
in each cell. In general, however, an insulating crystal 
is not a sum of isolated units, and there is an explicit 
dependence upon P mac due to transfer of electrons and 
derealization of the x hole between cells. 

There are other interesting conclusions which follow 
from the analysis of the periodic system in the presence 
of an electric field, e.g., the independent-particle kinetic 
energy Ti n d is also a function of the polarization . The 
generalized KS functional can be written Q| 

E[n; P mac ] = E ext [n; P mac ] 

+T in d[n; P mac ] + E H [n] + E xc [n; P mac ] , (4) 

where the first term is given in Eq. ([!]). The variational 
principle requires that E be stationary w.r.t. both n(r) 
and P mac , and in terms of the lattice-periodic single- 
particle orbitals ^ t of the KS approach, the Euler- 
Lagrange equations become M 



-i(V + ik) 2 + V eff (r) + i E eff • V k 



<t>i,k(r) 
£i,k(r) 0i,k(r) ■ 



(5) 



Here V eff (r) = v ext (r) + S(E H +E xc )/Sn(r) is an ordinary 
local potential and E e ff = E mac — j^dE xc /dP nmc is an 
effective electric field. Accordingly, the resulting gener- 
alized KS equations are formally the same as the usual 
density- only ones except that there is an additional term 
involving the derivative of the periodic orbitals with re- 
spect to the wave vector k (where k e [— ^; -1 in Id). 

The term dE xc /dP mac = —Q E xc in the generalized 
KS equations can be identified with an "x-c electric field" 
in analogy to the actual electric field term. It must 



be emphasized, however, that this is not a true electric 
field; it acts only on electrons and not on other charges; 
it may be longitudinal or transverse and is not related to 
any density by a Poisson-like equation. Since the changes 
are quadratic in E mac and thus also in P mac , it follows 
that in the case of a symmetric crystal E xc oc E mac for 
small E mac . Thus for zero field the system is described 
by the usual density-only KS equations, but E xc affects 
the dielectric properties. We do not attempt to estimate 
the size of the effect since the present Id model does not 
fix the full 3d form of the x hole needed to calculate the 
x energy. 

If one considers a system with lower symmetry, in gen- 
eral the ground state polarization P mac is non-zero, and 
there is no symmetry reason for E xc to vanish. Our cal- 
culations also provide a microscopic basis for such a field. 
If we consider aid example in which the crystal potential 
is non-mirror-symmetric, the shape of the x hole changes 
linearly with the field E mac . An example of the change 
in the hole per unit field is shown in Fig. |[ Since P mac 
varies linearly with E mac , it follows immediately that 
E xc 7^ 0, even for E mac = 0. Note that the sign of the 
change is relevant, and the positive direction for E mac in 
Fig. H is chosen to be in the same direction as P mac . The 
existence of a spontaneous E xc was discussed in a recent 
paper by GGG M, who described the counter- intuitive 
consequences if one attempts to incorporate it into the 
original density-only KS formulation. Vanderbilt |9| went 
on to find simple examples where straightforward appli- 
cation of the usual KS approach leads to "ultra-nonlocal 
dependence upon the charge density". In contrast, we 
have shown that the generalized KS framework provides 
a simple, direct description in terms of the dependence 
of the x-c hole upon the polarization as well as density. 

It is interesting to note that the present theory of 
polarized insulators is formally analogous to the well- 
known case of spin density functionals, in which a spin- 
polarized system is described by KS equations with a 
spin-dependent effective potential even if there is no such 
magnetic-field-like terms in the original HK system. It 
is necessary to allow for symmetry breaking terms in the 
KS functional in order to describe a broken symmetry 
state. A general KS formulation involving properties 
other than the density has been developed by Jansen jl5| , 
who showed that such a theory must in general involve 
operators in the effective Hamiltonian which couple to 
the desired additional variable(s). Jansen's analysis leads 
to equations similar to our Eq. (|^) in which the effective 
field in principle can be defined in terms of the exact 
many-body wave functions. 

In conclusion, we have provided a simple proof that 
in an insulating crystal, the dielectric properties cannot 
be described solely in terms of the periodic bulk density. 
This leads to a generalized KS approach needed for all 
materials with non-zero polarization, in which there is an 
effective "x-c electric field" E xc that can be represented 
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as a non-local operator. Our central result is that the x-c 
hole itself is modified in a fundamental way in a polar- 
ized medium. This was explicitly demonstrated for the 
x hole in a non-interacting system, and it follows from 
well-known expressions involving a coupling constant in- 
tegration H that the same conclusions apply to the x-c 
energy E xc in the interacting system. Our conclusions 
have no effect upon work done within approximations 
such as Local Density (LDA) and Generalized Gradient 
(GGA); instead, this work shows that a complete de- 
scription of dielectrics is outside the framework of such 
approximations. 
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FIG. 1. Results (in a.u.) for a model of a non-polar 
insulator with the lowest band occupied, i.e., 2 elec- 
trons/cell. The potential energy (see text) is given by 
Vo{x) = -V a cos{2nx/a), V a = 1/4, and V b = V a /5, with 
a = 5ao (ao = U 2 /(m e e 2 )). We show the central 8 cells 
in a supercell of length L — 80a. (L — 40a gives essen- 
tially identical results.) Top panel: V(x); the average slope 
is E mac = 2nVt/L. Middle panel: Change in electron density 
An(x) — n(x) — n(x)v b =o in the supercell, and in the fit- 
ted system (indistinguishable). Lower panel: Change in the 
density- weighted average x hole divided by |E mac | 2 for two 
different potentials which give the same periodic density. 



FIG. 2. Same as the lower panel in Fig. 1, but for a po- 
lar insulator with Vo(x) = —V a [cos(27r:r/a) — i sin(47rx/a)l , 
V a = 1/4, and V b = V4/20, plotted per unit field E mac . The 
spontaneous polarization is P^, ac = 1.321 x 10~ 2 , while the 
changes are SP^L = 0.86 x 10" 3 , and 5P%L = 0.52 x 10" 3 . 
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Figure 1: Gerardo Ortiz, Ivo Souza, and Richard M. Martin 




Figure 2: Gerardo Ortiz, Ivo Souza, and Richard M. Martin 



